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Abstract 

Let W be a correlated complex non-central Wishart matrix defined through 
W = X^X, where X is n x m (n > m) complex Gaussian with non-zero 
mean T and non-trivial covariance S. We derive exact expressions for the 
cumulative distribution functions (c.d.f.s) of the extreme eigenvalues (i.e., 
maximum and minimum) of W for some particular cases. These results are 
quite simple, involving rapidly converging infinite series, and apply for the 
practically important case where T has rank one. We also derive analogous 
results for a certain class of gamma- Wishart random matrices, for which 
Y H Y follows a matrix-variate gamma distribution. The eigenvalue distri- 
butions in this paper have various applications to wireless communication 
systems, and arise in other fields such as econometrics, statistical physics, 
and multivariate statistics. 
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1. Introduction 



Eigenvalue distributions of Wishart random matrices arise in many fie 



Prominent examples include wireless communication systems 



41 
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21 
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26 



ds. 



39 



53| . synthetic aperture radar (SAR) signal processing [32J, econo- 
metrics 51 1 ; statistical physics (4, 55], and multivariate statistical analysis 



48|. In many cases, the Wishart matrices of interest are com- 



plex [18], correlated, and non-central. Such matrices arise, for example, in 
multiple-input multiple-output (MIMO) communication channels character- 
ized by line-of-sight components (i.e., Rician fading) with spatial correlation 



amongst the antenna elements 



39] 



In this paper, a main focus is on the distributions of the extreme eigen- 
values (i.e., maximum and minimum) of Wishart matrices, which arise in 
many areas. For example, in the context of contemporary wireless communi- 
cation systems, the maximum eigenvalue distribution is instrumental to the 
analysis of MIMO multi-channel beamformin g sy stems 25] and the analysis 



of MIMO maximal ratio combining receivers [26|, |4l| , whereas the minimum 
eigenvalue distribution is important for the desi gn a nd analysis of adaptive 
MIMO multiplexing-diversity switc hing systems 22j, as well as the analysis 



of linear MIMO receiver structures 



44J. In the context of econometrics, the 



minimum eigenvalue of a non-central Wishart matrix is important for char- 
acterizin g th e weak instrument asymptotic distribution of the Cragg-Donald 



statistic 



51] . In statistical physics, information pertaining to the nature of 
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entanglement of a random pure quantum state can be obtained from the 



two extreme eigenvalue densities of Wishart matrices 



Moreover, the 



maximal and minimal height distributions of iV non-intersecting fluctuating 
interfaces at the thermal equilibrium and with a certain external potential are 
also related to the extreme eigenvalues of a Wishart matrix [43|. As a final 
example, in SAR signal processing, the probability density of the maximum 
eigenvalue a Wishart matrix is an important parameter for target detection 



and analysis 



22]. 



We focus primarily on correlated complex non-central Wishart matrices, 
as well as another important and closely related class of random matrices, 
which we refer to gamma- Wishart. Such matrices arise in the context of 
MIMO land mobile satellite (LMS) communication systems jl|, and corre- 
spond to non-central Wishart matrices with a random non-centrality ma- 
trix having a distribution which is intimately related to the matrix-variate 
gamma. As discussed in [1], the eigenvalues of gamma- Wishart random ma- 
trices are important for the design and analysis of MIMO LMS systems; 
for example, the maximum eigenvalue density determines the performance 
of beamforming transmission techniques, whereas the minimum eigenvalue 
density is closely related to the performance of linear reception techniques. 

Recently, the marginal eigenvalue distributions of random matrices have 



received much attention; for surveys, see 



y 



54( . For the extreme eigen- 



values, distributional results are now available for correlated central, uncor- 
rected centraL and uncorrelated noncentral complex Wishart matrices (see, 



for example, 



131, 



16 
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41 
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m F» less is 

known for gamma- Wishart matrices, other than the results in [1] , which deal 
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exclusively with uncorrelated matrices. In the majority of cases, the stan- 
dard approach has been to integrate the respective joint eigenvalue densities 
over suitably chosen multi-dimensional regions. For the more general class of 
complex non-central Wishart and gamma- Wishart matrices with non-trivial 
correlation however, there appears to be no tractable existing results. For 
these matrices, as we will show, the joint eigenvalue densities are extremely 
complicated, and it seems that this direct approach cannot be easily under- 
taken to yield meaningful results. 



In this paper, by emplo ying 



n 



[fly, 



an alternative derivation technique (also con- 



sidered in [10|, LLU |30|, [37|, |42|, |46() which allows us to deal with the joint 



matrix- variate density rather than the density of the eigenvalues, we derive 
new exact expressions for the cumulative distribution functions (c.d.f.s) of 
the minimum and maximum eigenvalues of correlated complex non-central 
Wishart and correlated gamma- Wishart random matrices. In both cases, 
whilst a general theory which accounts for all matrix dimensions and distri- 
butional parameters appears intractable, we are able to derive solutions for 
various important scenarios. Specifically, for correlated non-central Wishart 
matrices, we derive expressions for the minimum eigenvalue c.d.f.s when the 
matrix dimensionality and the number of degrees of freedom are equal. We 
also derive results for some specific scenarios for which they are not equal, 
and present some analogous results for the maximum eigenvalue c.d.f. For 
tractability, we focus on matrices with rank-one non-centrality parameter, 
which is practical for various applications; most notably, MIMO communi- 
cation systems with a direct line-of-sight path between the transmitter and 
receiver. Given the overwhelming complexity of the underlying joint eigen- 
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value distribution, these extreme eigenvalue c.d.f. expressions are remarkably 
simple, involving infinite series with fast convergence, and they can be easily 
and efficiently computed. 

For the case of gamma- Wishart matrices, we focus on scenarios for which 
the underlying matrix-variate gamma has an integer parameter. The implica- 
tions of this assumption from a telecommunications engineering perspective 
are discussed in jl|. As for the non-central Wishart case, we derive exact ex- 
pressions for the minimum and maximum eigenvalue distributions for certain 
gamma-Wishart particularizations. 

Whilst previous expressions pertaining to the non-central Wishart case 
have been reported in [if} 0, 46]; those are very complicated, involving 
either infinite series' with inner summations over partitions with each term 
involving invariant zonal polynomials (c.f. Section [2]), or infinite series with 
special functions of matrix arguments jll[ 3^]. As such, those previous results 
have limited utility from a numerical computation perspective. 



2. Preliminaries and New Matrix Integrals 

2. 1 . Preliminaries 

In this section, we provide some preliminary results and definitions in 
random matrix theory which will be useful in the subsequent derivations. The 
following notation is used throughout the paper. Matrices are represented as 
uppercase bold-face, and vectors by lowercase bold-face. The superscript {-) H 
indicates the Hermitian-transpose. I p denotes apxp identity matrix. We use 
| • | to represent the determinant of a square matrix, tr(-) to represent trace, 
and etr(-) stands for exp (tr(-)). The set of complex Hermitian mxm matrices 
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are denoted by "H m and the set of Hermitian positive definite matrices are 
denoted as For A, B G 7i m , A > is used to indicate the positive 

definiteness, and A > B denotes A — B G A > is used to indicate 
non-negativeness. Aj^ represents the j, kth element of matrix A. \x] is 
the ceiling function, defined as \x] = min {n G Z|n > x}. Finally, the kth 
derivative of function f(y) is represented as f^ k \y) for all k G Z + , and with 

f m {y)-.= f{y). 

Definition 1. The generalized hypergeometric function of one matrix argu- 
ment can be defined aj^ 

[ai] K [a2] K ---M«C«( Y ) 

k=0 K 



P F ? (oil a 2 , . . . , a p , 6 l5 b 2 , . . . , b q ; Y) - ^ ^ _ _ _ _ 



where Y G "H m; [a] re = Y\_( a — i + -Ofcj; K = ^2, • • • , &m) partition 

3=1 

of k such that ki > k 2 > . . . > k m > and YllLi ^» = ^, an< ^ ( a )fc = 
a(a + 1) • • • (a + — 1). yMso, £/ie complex zonal polynomial C K {Y) is defined 



in 



MJ 



Remark 1. Note that the infinite zonal polynomial expansion given in (QJj 
reduces to a finite series if at least one of the ajS is a negative integer. As 
such, when N G Z + we have 

"' A - : N] K [a 2 } K ---[a p ] K C K (Y) 



-N, a 2 , . . . , a p ; b u b 2 , . . . , 6 g ; Y) = ^2 K ~\f 

1 — n L 1 



(2) 

where Y1 K denotes the summation over all partitions k = k 2 , . . . , k m ) of 
k with k\ < N. 



1 The convergence of the infinite zonal series is discussed in 



For more properties of zonal polynomials, see [5|, |24j, |52|. 



Definition 2. (Non-Central Wishart Distribution) Let X be an n x m 

(n > m) random matrix distributed as CAI ' n<m (T,I n <g> S) ; where £ £ 
and T £ C nxm . Then W = X H X £ H+, has a complex non-central Wishart 



distribution W m (n, S, 0) with density function l29l ] 



P f r (—ft) |w|"~ m ~ 
/w (W) = ^ ; etr (-S-W) F X (n; GS^W) (3) 

where = E -1 T H T zs t/ie non-centrality parameter and T m (-) represents 
the complex multivariate gamma function defined as 

m 

~ * A m(m — 1) i — r , * 

T m (n) = n^- Y[T(n-j + l) 

3=1 

with T(-) denoting the classical gamma function. 

Definition 3. (Matrix Variate Gamma Distribution) Let a > m and 

fi £ The random matrix M £ /ias a matrix-variate complex gamma 
distribution T m (a, f2) i/ zfo density is 13a , Def. 6.3]. 

Definition 4. (Gamma- Wishart Distribution) Let us construct annxm 
matrix X such that 

X = X + X (4) 

where X ~ CA/ ' n<m (0, I n ® S) and X X ~ T m (a, ft) are independent. Then 
V = X H X £ follows a gamma- Wishart distribution rW m (n, a, X, f2) 
gwen fry ^/ 

, . etr (-E _1 V) |V| n - m |fJ| a ~ / w , . _i 1 \ 

/V(V) = fluEi-l^ + np lF ' ( ai n; E < E + n > s v ) ' 

(5) 

iVote i/iai /or a = n, (TJ]] reduces to W m (n, S + f2 _1 ) . 

7 



In addition to zonal polynomials, non-central distributional problems in 
multivariate statistics commonly give rise to other classes of invariant poly- 
nomials 



The next lemma presents the joint eigenvalue distributions o 
Wishart matrix, in terms of invariant polynomials defined in 11 



gamma- 



12 



46]. 



the correlated 



4& Eq. 5.4] and 



The proof of this lemma follows similar steps to the proo 
non-central Wishart joint eigenvalue density, g\ (A), in 
thus omitted. 

Lemma 1. The joint density of the ordered eigenvalues Ai > A2 > ■ ■ • > 
A m > 0, of the matrix V in (TJJ) is given by 



0a(A) 



r m (n)r m (m)|£Mft + £ 



k=l 



k<l 



k,s=0 



[<cj CT (-s- 1 , s^ 1 (n + s- 1 )- 1 s- 1 ) (A, A) 



k\s\[n] a C^l m ) 



(6) 



where A is a diagonal matrix containing the eigenvalues 0/ V along the main 
diagonal. 



The following technical lemma is proved in |Appendix A 



Lemma 2. Let Xi,x% be the two distinct eigenvalues o/X G H2 ■ Then, for 
all n G Z + , 



IV1 

^ (-l) i 4*e l ri |X| i tr n - 1 - 2i (X) 



(7) 



i=0 



where ef denotes the ith elementary symmetric function of the parameters 

~n-2 



S n := Wf-Vcos 2 ( — ) cos 2 



n 



n 



(8) 



2.2. New Matrix Integrals 

Here we present some new matrix integral results which will be impor- 
tant in the derivations of the extreme eigenvalue distributions, given in the 
following sections. 

Lemma 3. Let A G H 2 and B G "H 2 with B > 0. Also, define Xi(y) and 
x 2 (y) as the eigenvalues of A + By. Then, Vp G Zq and 3?(a) > 1, 



f 

Jo 



X| a ~ 2 etr (AX) tr p (BX) rfX = ^ 



f 2 (a)f 2 (2) 



r 2 (a + 2) 



; A,B 



.(0) 



(9) 



where 0a?b,o(O) is calculated recursively via 

v 



«B, a (°) = j"^) ( A 2ka(0) " E ^ 



with initial condition 



0i°ka(O) = 0A,B,a(O) 



A A ,B,a(0) 

x 1 (0)-x 2 (0) 



(10) 



;n) 



Here, 



Aa,b,o(j/) = a;i(j/)iFi (a; a + 2; Zi(y)) i-Fi (a - 1; a + 1; x 2 (y)) 

- x 2 (y) 1 F 1 (a; a + 2; x 2 (y)) iF : (a - 1; a + 1; Xl (y)) (12) 



and 



with 



x[ j) (0) 



h { i] B (0)=x?(0)-x«\0), 

' xi(0)tr(B) - |A|tr (BA^ 1 ) 
x 1 (0)-x 2 (0) 
2(4 1) (0)4 1) (0)-|B|) 

x 1 (0)-x 2 {0) 

EU(D^" fc) (Q)4 fc) (o) 

x 1 (0)-x 2 (0) 
9 



(13) 



if j = l 

if j = 2 
if j>3, 



(14) 



lAltrfBA- 1 ) -z 2 (0)tr(B) . = 

„(j)/~N I _ /nN _ /rA 1 J — 



a#'(0) 



si(0) - x 2 (0) 
-x[ j) (0) 



(15) 



if j>2. 



Proof. See Appendix B 



Lemma 4. Let A G "H+ and let R G "H m twt/i unit rank. Then, for t G Zq 
and 9ft(a) > m — 1, 



etr (—AX) tr (X) |X| a_m tr* (RX) dX 



)*f m (a)tr* (RA- 1 )^- t 



tr (R (A- 1 ) 2 
tr(RA- 1 ) / +atr(A ' 1) 



(16) 



Proof. See Appendix C 



When the matrices are of size 2 x 2, we can obtain the following general 
result: 

Lemma 5. Let A G H% and let R G H2 with unit rank. Then, for p, t G Zq 
and 3?(a) > 1, 



[ etr (-AX) tr p (X) |X| a ~ 2 tr* (RX) dX 
r (a) t f 2 (a) (-l) fc (D 



£ ^^tr*-*(RA->*(R)C^ 



tr (A) 



2a/|A| 



(17) 



where C^(-) denotes an ultraspherical (Gegenbauer) polynomial. 



Proof. See Appendix D 
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Lemma 6. Let A G 71% and let R(> 0) G T-L% with unit rank. Then, for 
t G Z+ ; 



I etr(-AX)tr i (RX)C ljli0 (X)rfX 



f 3 (4)|A|- 4 ((4) t tr* (RA- 1 ) tr(A) + t(4),_ 1 tr'- 1 (RA" 1 ) tr(R)) . 

(18) 



Proof. See Appendix E 



Lemma 7. Let A, B G T-L^ ■ Then, for p,t G Zq and 3?(a) > 1 ; we have 
[ etr (-AX) tr p (BX) tr* (X) |X| a ~ 2 dX 

= P m|A| r 2 (a) £ (tl + 1)!(t _ tl)! L s ™ 



i=0 



(19) 



where 



min(p, £ tl>i ) 



fc=0 ^ ' 



p — k p — k 



e *i.«- fc (A)tr fc (B) |A|- £t i- V |B 
£tl+a / tr (A^B 



x C 



p—k 



2y^B\ J ' 



et^i = 2ti — * — 2i, Et 1 =t\—i, and r = (ti,t — *i) is a partition oft such that 
[§] < *i < *■ Moreover, e[ denotes the ith elementary symmetric function 
of the parameters 



S T := < cos 



7T 



2*i - * + 1 



, cos 



2tt 



2*! - * + 1 



. . . , COS 



2ti — t — 1 



7T 



2*i - * + 1 



(20) 
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Proof. See Appendix F 



Armed with the new results in this section, we are now in a position to 
derive the extreme eigenvalue distributions of both correlated complex non- 
central Wishart and gamma- Wishart matrices. These key results are the 
focus of the following two sections. 

3. New Minimum Eigenvalue Distributions 

In this section, we consider the minimum eigenvalue distribution. To 
evaluate this, the most direct approach is to integrate the joint eigenvalue 
probability density function (p.d.f.) as follows: 

F min (x) = 1 - P(Ai > ■■■> X m > x) 

= 1-1 g(\)dX 1 ---dX m (21) 

where V = {x < X m < ■ ■ ■ < Ai} and g(A) G {(7a(A), g^(A)}. This direct 
approach, however, is difficult for two main reasons: (i) due to the presence 
of the invariant polynomials in the joint eigenvalue densities, and (ii) due the 
unbounded upper limit of the integrals which makes term-by-term integration 
intractable. To circumvent these complexities, in the following we adopt an 
alternative derivation approach based on integrating directly over the matrix- 
variate distribution itself, rather than the distribution of the eigenvalues. 

To highlight the approach, consider Y G with minimum eigenvalue 
Amin(Y) having c.d.f. 

F min (x) = P (A mm (Y) < x) = 1 - P (A mm (Y) > x) . (22) 
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The key idea is to invoke the obvious relation^ 

P (X min (Y) > x) = P (Y > xl m ) (23) 

which allows one to deal purely with the distribution of Y, rather than the 
distribution of its eigenvalues. 

3. 1 . Correlated Non- Central Wishart Matrices 

For the non-central Wishart scenario, we deal with the matrix W with 
joint density given in (J3J). Thus, with f l2"3"j) . we have 

P (WW) > z) = / /w (W) dW 

JW>xI m 



exp (—77) 
r m (^) |S| n Jw-xi m eH 



[ |W| n - m etr (-S^W) 

x Fi (n; 0S _1 W) dW (24) 



where 77 = tr(0). Applying the change of variables W = x (I m + Y) with 
dW = x m2 dY yields 



j,mn , 



exp (— 77) etr (— x£ x ) f 

P (A mm (W) > x) = V ; " \ n } - / |I m + Y| 

Y m (n) £ ivew+ 



x 



etr (-xS _1 Y) Fi (n; (I m + Y)) cTY 



It is convenient to now expand the hypergeometric function with its equiva- 
lent zonal polynomial series expansion ([I]) to give 

P(x minm > X ) = ^^7\t, ( r s " 1} ee 1 



x / |I m + Yr" m etr(-sS- 1 Y)a(x©S _1 (I m + Y))dY (25) 
This relation has also been employed previously in 

, 0,0, 0,0,0- 
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where k = K m ) is a partition of k into not more than m parts such 

that «i >....> K m > and Ki = k. 

Observing that @E 1 is Hermitian non-negative definite with rank one, 
it can be represented via its eigen decomposition as 

0£ 1 = ficta 11 (26) 

where a e C mxl and ol h ol = 1. Recalling that zonal polynomials de- 
pend only on the eigenvalues of their matrix arguments, and noting that 
©X 1 (I m + Y) is also rank one, we can write 025p with the aid of f[2"6"j) as 

P(A-(W) > ,) = I "" eXP '7\tt" XS "' ) EE ' 



X 



/ |I m + Y| n " m etr (-xS^Y) C K (^a ff (I m + Y) a) dY. (27) 



Applying the complex analogue of [42j, Corollary 7.2.4], since a. H (I m + Y) ck 
is rank one, then it follows that C K [xjicx. 11 (I m + Y) a) = for all partitions 
k having more than one non-zero part. Hence 

C K (x^ H (I m + Y) a) = (x^) k (f) ti* (aa H Y) (28) 

t=o ^ ' 

and (|2"7|) can be written as 

r m (n) |S| W 

(29) 

where 

Q^ n (x)= / |I m + Yr- m etr(-xS^Y)tr' (q^Y)^. (30) 

14 



Unfortunately, it appears that this integral is not solvable in closed form for 
arbitrary values of m and n. However, as we now show, it can be solved in 
closed-form for various important configurations, thus yielding exact expres- 
sions for the minimum eigenvalue distributions. These results are presented 
in three key theorems. In each of these, we recall the notation 

p = tr (eS -1 ) , 77 = tr(0). (31) 

The theorem below gives the exact minimum eigenvalue distribution for 
"square" Wishart matrices: 

Theorem 1. Let X ~ CN m>m (T,I m £g> S) ; where T e C mxm has rank one, 
and W = X^X. Then the c.d.f. of A m j n (W) is given by 

F mm {x) = 1 - exp (- V ) etr (-xS" 1 ) V j^^iFx (m; m + k; V ) . (32) 

k\[m)k 

Proof. Substituting m = n into fl29|) and fl30l) yields 



(33) 



where 



Qljx) = [ etr (-xS^Y) C T (aa fl Y) dY. (34) 
This matrix integral can be solved using 38|, Eq. 6.1.20] to give 

Ql m = fm(w)(m > |sr a f -) = fm(m)(m 2 )t|s|m f -5-1 * (35) 

m ' mV ' x m \/ix ) x m \xfi ) v ' 
where we have applied (126]) to arrive at the argument of the zonal polynomial. 

Substituting ( 1351) into ( |33i) with some manipulation yields 

P (A mm (W) > x) = exp (- V ) etr (-xS" 1 ) £ £ f ' . 



(36) 
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To obtain a power series in x, we re-sum the infinite series as follows 

EKe(') (™)« (-3-V = E to m + *i 5) ■ (37) 

Finally, using (13"?]) in (I55]l with ((22]) gives the result in (pSj). □ 

Remark 2. v4n alternative expression for the c.d.f. can be obtained by ob- 
serving the fact that 

= $3 (m, m, rj, xfj,) (38) 

where $>s(a,b,x,y) is the confluent hypergeometric function of two variables 
la . Eq. 5.7.1.23 }. Thus, we can write the minimum eigenvalue c.d.f. as 

F m in{x) = 1 - exp (-rj) etr (-xE -1 ) $ 3 (m, m, r], xfx) . (39) 

The theorem below gives the exact minimum eigenvalue distribution for 
2x2 Wishart matrices with arbitrary degrees of freedom: 

Theorem 2. Let X ~ CM n ,2 (T, I n <g> £), where T G C nx2 Aas rank one, 
and W = X H X. JTien i/ie ad./, o/ A m j n (W) is given by 

etr (-xE- 1 ) ~ (x/x) fc ^(k\( V ^ t 



(40) 



w/iere 

n-2 i min(j',t) , n \ / -\ /jA / \ ' 



») = EE E (-1)' (" 7 2 ) (;) (««) j . (-.,) (£) 



and = i — j + 2. 
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Proof. We begin by substituting m = 2 into (1211 and (150j) to yield 

P (W W) > x) = ^^x-etr (-xS-) GU*)- 
r 2 (n)|S| ^/c!(nj fe ^VV 

Now we use the determinant expansion 

|I 2 + Y|- 2 = g £ f " T ^ H ^ (Y) |Yf- (41) 



to write Q% „ (x) as 



i=0 i=l 



Gl»(*) = EE( n /) ( •) / YeW9+ trJ oo i vr-w (-xs-y) 

tr* (aa^Y) dY. (42) 



X 



Lemma O can be used to solve the above integral in closed form and subse- 
quent use of (1221) followed by some algebraic manipulations gives (j4"Dl) . □ 

Although the c.d.f. result in Theorem [2] is seemingly complicated, it can 
be evaluated numerically for any value of n. Moreover, for specific values of n 
it often gives simplified solutions. Some examples are shown in the following 
corollaries. 

Corollary 1. Let X ~ CAf 3)2 (Y, I 3 ® S) ; where T G C 3x2 has rank one, 
and W = X^X. Then the c.d.f. of \ min (yV) is given by 

F mm {x) = 1 - exp (-7?) etr (-x^ 1 ) V g^- ^(fc, 77, x) (43) 

w/iere 

J* 3)2 (A;, 77, x) = ft(x) x Fi (3; 3 + k; rf) + ^2(^)1^1 (2; 3 + k,rj) , 

Qi(x) = 1 + ( tr (XT 1 ) ) x, and ^(x) = -x + — — . 

\ 7) J rj 2\2j\ 
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Remark 3. An alternative expression for the above c.d.f. can be written 
based on the confluent hypergeometric function of two arguments as 

F m in(x) = 1 - exp (-77) etr (-xE -1 ) (ai(x)$ 3 (3, 3, 77, a;^) 

+ a 2 (x)$ 3 (2, 3,77,^)) ■ 

Corollary 2. Let X ~ CM ^2 (Y , I 4 ® S), 7t;nere T e C 4x2 /ias ranA; one, 
and W = X^X. Then the c.d.f. of A min (W) is given by 

F min (x) = 1- exp (-77) etr (-xE^ 1 ) ^n - * 774 - 2 ^' ?7 ' ^ ^ 

fc=0 ^ fe 



where 



«, a;) = ^i(aOiFi (4; 4 + fc; 77) + i^(x)iFi (3; 4 + fc; 77) 

+ z/ 3 (x) 1 F 1 (2;4 + fc;77), 



/ \ 1 a i 2 

1/1(2; J =l + aix + — x , 



u 2 (x) =^x + Q + y + jjtr (S x ) ai - a 2 ^) x 2 + 



Jb 



= (^- ktr - £ + ^ + b ^ 



a? 



3|Sf ' 

2 



6 



ux 3 x 4 

+ + 



3t7|S| I2|£| 2 ' 
ai = tr (S- 1 ) - * and a 2 = tr 2 (S" 1 ) - ^ - ji. 

Remark 4. An alternative expression for the above c.d.f. can be written as 

F min (x) = 1- exp (-77) etr (-xE" 1 ) (i/i(x)$ 3 (4, 4, 77, x/i) 

+ z/ 2 (x)$ 3 (3, 4, 77, x/i) + z/ 3 (x)$ 3 (2, 4, 77, x/i)) . 
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The theorem below gives the exact minimum eigenvalue distribution for 
3x3 Wishart matrices with 4 degrees of freedom: 

Theorem 3. Let X ~ CjV 4 , 3 (T, I 4 ® S), where T G C 4x3 has rank one, 
and W = X H X. T/ien the c.d.f. of A m j n (W) zs given by 

F mm {x) = 1 - exp (-77) etr (-X5T 1 ) £ ^T 7 " 4 ' 3 ^' ^ x) (45) 

fc=o fc 'W* 

J" 4)3 (A;, 77, x) = pi(x) 1 F 1 (4; 4 + fc; 77) + p 2 (^)i^i (3; 4 + k; rj) , 

^ _ , tr(0S) 2 



pi{x) = 1+ ftr (IT 1 ) + 
p 2 (x) = -x + — ^- ftr (£) - tr (0S) x 2 + 



2r7|S| 



77 2|S| V 77/ 6|S|' 

Proof. We can write ( 1291) and ( 1301) in the case of m = 3 and n = 4 as 

^ '■■ \ r, 

'3,4 



P ( W W) > „ = g^etr (-^) g g £ (J QUx) 



(46) 



Following the identity 

|I 3 + Y| = 1 + tr(Y) + |Y| + Ci,i, (Y), (47) 
we can write <2 34 (x) as 

Q\ A (x) = [ etr (-x5T x Y) tr* (ctct H Y) dY 

+ [ etr (-xS _1 Y) tr(Y)tr* (q^Y) rfY 

+ / etr (-xE^Y) |Y|tr* (aa ff Y) dY 

+ / etr (-xS _1 Y) C lil)0 (Y)tr* (aa ff Y) dY. (48) 
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These matrix integrals can be solved with the aid of 38|, Eq. 6.1.20], Lemma 
HI and Lemma |6] to yield 



i ? . i ■■■ i . , i 11 / I ) v 



|S| 4 f 3 (4) 


(-) 


X 12 


\XflJ 



QUx) = ' ' 12 av ' ± (Pi(^)t + P2(x)(3) t ) (49) 



where we have used the relations t(3) t = 3(4) t — 3(3) t and £(4) t _i = (4) t — (3) t . 
Substituting (j4"9l into f|46l) . we obtain 



P(A mm (W) >x) = exp(-^)etr(-xS- 1 ) L (4)/^-) ' 

\ fc=o ' k 4=0 ^ ' V /V 

Finally, we re-sum the infinite series as power series in x and use ( 1221) to 
arrive at the result in (USD. □ 



Remark 5. v4n alternative form of the c.d.f. above can be written as 

F min {x) = 1 -exp(-r)) etr (-xS _1 ) (pi(z)$ 3 (4, 4, rj, x/i) 

+p 2 (x)® 3 (3,4,r],xiJ,)) . 

We now present some simulation results to verify the validity of our new 
minimum eigenvalue distributions. We construct the covariance XI matrix 
with (j, k)th element 

^ = exp^-^(j-A;) 2 ^, l<j,k<m (50) 

and the mean matrix T as 

T = a^b (51) 
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where 



a = [1 exp (2m cos 9) exp (Am cos 9) . . . exp (2(n — l)iir cos 9)} 
b = [1 exp (2m cos 9) exp (Am cos 9) . . . exp (2(m — \)m cos 9)] 

with # = 7i /A and i = V - 1- Note that these particular constructions for 
the covariance and mean matrices are employed since they are reasonable for 

nn 

modeling practical correlated Rician MIMO channels |3|, |39[ . 

Fig. [1] compares our analytical results with simulated data. The analyti- 
cal curves for the cases m = n were calculated based on Theorem [TJ while for 
the cases m = 2 and m = 3, they were calculated based on Theorems [2] and 
[3] respectively. The accuracy of our results is clearly evident from the fig- 
ure. Note that in evaluating these analytical curves, the infinite summations 
in ( 132|) . ( H0|) . and (H5i) were truncated to a maximum of 20 terms; thereby 
demonstrating a fast convergence rate for each series. 

3.2. Gamma- Wishart Matrices 

We now turn to the analysis of the minimum eigenvalue distribution of 
gamma- Wishart random matrices. In this case, we deal with the matrix V 
with joint density given in ([3]). Thus, with ( 12"3"|) . we have 

P(Xmvn(V) >x) = K m>n [ |V| n ~ m etr(-S- 1 V) 1 F 1 (a;r2;SV)dV 
where S = S" 1 (S" 1 + Q)' 1 S" 1 and /C m , n = fm(n)| J"i'°- 1+nr • Applying 



the change of variables V = x (I m + Y) and using the Kummer relation 23] 
1F1 (a- n- xS (l m + Y)) = etr (xS (l m + Y)) X F X (n - a; n; -xS (I m + Y)) 
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0.5 1 1.5 2 

X 



(a) n = m 




x 

(b) to = 2,3 

Figure 1: Comparison of the analytical minimum eigenvalue c.d.f.s with simulated data 
points for correlated non-central Wishart matrices of various dimensions. 
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yields 

P (A mm (V) >x) = JC m , n x mn etr (-xQ) [ \I m + Y^etr (-xQY) 

x 1 F 1 (n-a;n;-xS(I m + Y))dY (52) 

where Q = £ 1 - S. 

This integral seems intractable for arbitrary values of m, n, and a. How- 
ever, as we now show, it can be solved in closed form solutions for some 
important configurations, thus yielding new exact expressions for the mini- 
mum eigenvalue distributions. 

The theorem below gives the exact minimum eigenvalue distribution for 
2x2 gamma- Wishart matrices with arbitrary degrees of freedom (i.e., arbi- 
trary n). 

Theorem 4. Let V ~ TW 2 {n, a, S, ft), with a e Z+ such that a > n > 2. 
Then the c.d.f. of A min (V) is given by 

r2fc 1 -fe-l-| 

2(a— n) min(fc,a— n) \ 2 I 

F mm (x) = l-JC 2 , n x 2n etr(-xQ) £ £ ]T dfl kul (x)x k 

k=0 fcl= |"fc-| «=0 

(53) 

fcl _ (a-n)!(a-n + l)!(2fci-A; + l) 

1 ~ (a - n - fci)!(a - n + 1 + fci - k)\ (h + 1)! (k - - l) fc _ fel 

= (-l)^ef |S| fe - fel+ '. 



,^,Ak,A tr^^(S) ' j! 
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with 

tl =[|] i=0 

w /i ere 

min(p,e tll i) 



g=0 \ ? / 

f tr^S) 



« = (fci, — fci) zs a partition of k such that |~|] < ki < min(fc, (a — n)), 



r = (tx,t — ti) is a partition oft such that | < t\ < t, Ujj = j — t + 2 and 
Vki,i = n + l + k-k!-2. 

Proof. Particularizing fl52|) tom = 2, a>n>2 and a € Z + , and applying 
the zonal polynomial expansion (J2J) yields 

2(a-n). r , _ 

P(A mm (V) > x) = /C 2 , n :r 2 "etr(-xQ) V V [ ,, (-x) fc 

x / |I 2 + Yr- 2 etr(-xQY)a(S(I 2 + Y))rfY 

(54) 

where k = (fci, fc 2 ) is a partition of k into not more than two parts such that 
+ = k and fci > A; 2 > 0, V&i G {0, 1, . . . , a— n}. Note that the series over 
k is finite (truncated at k = 2(a — n)) due to the negative sign of the gener- 
alized complex hypergeometric coefficient. Careful inspection reveals that k 
can be written as k = (ki, k — ki), where [~~] < k\ < min (k, (a — n)). This 
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fact, along with the alternative representation of complex zonal polynomial 



given m 
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52j, and Lemma [21 



^( s ^ + Y ^ ( rfi)K^)! |S(l2 + Y)|fc " fel 



x Yl (-l) Z 4'ef|S(I 2 + Y)|V fe i>' (S(I 2 + Y)) (55) 

1=0 

gives (after some manipulations) 

2(a— n) min(fc,(a— n)) 

P(\min(V)>x)=JC 2 , n x 2n etr(-xQ) ^ £ d* 1 

x dfl kul (x)x k (56) 

z=o 

where 



[ etr(-xQY)|I 2 + Y| I/fc i' i tr £fc i. i (S(I 2 + Y))dY. (57) 
Using |I 2 + Y| = 1 + tr(Y) + |Y| and the binomial theorem yields 

™-%%<y){y)i&z%*±®w*^ (58) 

where 

J tpj = 1 , / etr(-sQY)tr p (SY)tr*(Y)|YP'-*dY. (59) 

Finally, solving the remaining integral using Lemma [7] and recalling ( l22l 

concludes the proof. □ 



Note that the minimum eigenvalue c.d.f. result given in ( |53|) can be easily 
computed numerically, since it contains only finite summations. Moreover, 
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for specific values of n and a, it leads to simplified solutions, as shown in the 
following corollary. 

Corollary 3. Let V ~ rW 2 (2, 3, S, SI) . Then the c.d.f. of A mm (V) is given 
by 

F mm {x) = 1- -Jp — etr(-xQ) (|l 2 + fT^ 1 1 + 

+ ^ + tr(Q- 1 )|S|)x+i|x 2 ). (60) 

The theorem below gives the exact minimum eigenvalue distribution for 
3x3 gamma- Wishart matrices with 3 degrees of freedom. 

Theorem 5. LetV ~ rW 3 (3, 4, S, SI). Then the c.d.f. of X min (V) is given 
by 

^( g ) =1 - l "s-i + n? ) ( |Is + frls " 1 + tr(F) i + tr(G)|s| y + |s| £) 

(61) 

where 

F = 2S-3CT 1 S-3|S|CT 1 + 6|S||Q|~ 1 Q + 3 |I 3 + S\ Q 1 (I 3 + S) _1 S (62) 
and 

G = S~ 1 + 3CT 1 . (63) 
Proof. In this case ( 152]) becomes 

P (A min (V) > x) = /C 3 , 3 x 9 etr (-xQ) / etr (-xQY) 

JYeHj 

x 1F1 (-1; 3; -xS (I3 + S)) tZY (64) 
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which upon applying the zonal polynomial expansion for the hypergeometric 
function fl2]) yields 

P (A mm (V) >x) = /C 3 , 3 x 9 etr (-xQ) V V / etr (-xQY) 

xC K (S{I 3 + Y))dY (65) 

where k = (^1,^2,^3) is a partition of k. It is not difficult to see that the 
admissible partitions corresponding to the integers 0, 1, 2, and 3 are (0, 0, 0), 
(1, 0, 0), (1, 1, 0), and (1, 1, 1) respectively. Thus, we can write (|65|) as 

P (A min (V) >x) = /C 3 , 3 x 9 etr (-zQ) f / etr (-xQY) dY 

+ ^ [ etr {-xQY) C 1)0 , (S (I 3 + Y)) dY 
+ I etr (-xQY) C M , (S (I 3 + Y)) dY 

[ etr (-xQY) C XX x (S (I3 + Y)) dY ) . (66) 



x 3 



Moreover, we have 

C 1A0 (S(I 3 + Y))=tr(S(I 3 + Y)) 

C 1 , 1 , 1 (S(I 3 + Y)) = |S||I 3 + Y|. (67) 

Utilizing (jJTj) we can express 

Ci.1,0 (S (I 3 + Y)) = |I 3 + (S (I 3 + Y))| - 1 - tr (S (I 3 + Y)) - |S (I 3 + Y) |. 

(68) 
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Now, using ( 157)) and (155)) in (155)) yields 

q , ^ (( x 2 tr(S)x 2 tr(S)a; 
P (A mi „(V) > ar) = /C 3 , 3 x 9 etr (-xQ) 1 - — - + 



6 6 



where 



and 



x / etr(-xQY)dY 

+ (f - t) f etr ( _X Q Y ) Ciao (SY) dY 
+i|l ( x 3 _ ^ ^ (a .) + )l3 + S | ^£ 2 (x)) (69) 

Gi(x)= f etr(-xQY) |I 3 + Y| rfY (70) 
./Ye«+ 

&(ac)= / etr(-a;QY) |l 3 + (I 3 + S)- 1 SY|rfY. (71) 



The first and second integrals in (169)) can be evaluated using |38|, Eq. 6.1.20], 
thus we concentrate on the evaluation of Q\{x) and Q2(x) . We provide a 
detailed solution for the integral Qiix) only, since both ( 170)) and ( 17T1) share 
a common structure. 

Using the relation |l 3 + (I 3 + S)" 1 SY| = y % (-1; - (I 3 + S)' 1 SY) in 
(171]) yields 

02(x) = [ etr (-xQY) x % (-1; - (I 3 + S)" 1 SY) dY. (72) 



This integral can be solved using 46j, Eq. 3.20] as 



02(x) = f 3 (3)|Q|- 3 x- 9 2 F (-1, 3; -x^Q- 1 (I 3 + S)" 1 S) 



3 ^-Qfc 



f 3 (3)|Qr 3 *- 9 £ -^rJ^-^lC, (Q- 1 (I 3 + S)- 1 S) . 

fc=0 ' 
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Since the valid partitions corresponding to the summation index k = 0,1,2 
and 3 are respectively (0, 0, 0), (1, 0, 0), (1, 1, 0) and (1, 1, 1), we can use equa- 
tions analogous to (1671) to obtain 

Q 2 (x) = f 3 (3)|Q|- 3 x- 9 (1 + 3x-hi (Q- 1 (I 3 + S)- 1 S) 

+ 6a;- 2 C lil ,o(Q- 1 (l3 + S)- 1 S) 
+ ear^lQl" 1 |Is + SI -1 |S|) . (73) 

Following similar arguments, we can obtain 

g x {x) = f 3 (3)|Q|- 3 x- 9 (1 + 3x-Hr (CT 1 ) + Qx- 2 C 1X0 (CT 1 ) + ex^lQ)- 1 ) . 

(74) 

Finally, using (1731 , (1741 , and (IE. II) in (|69|) , recalling ( 122]) , and applying some 
lengthy algebraic manipulations, we arrive at the result in (1611) . □ 

Fig. [2] compares our analytical results with simulated data. The analyti- 
cal curves for the cases m = 2 and m = 3 were computed based on Theorems 
|4] and [5] respectively. Here we have used the same £ as defined in (150]) , 
whereas f2 is constructed with the following j, kth element: 

/ 147^3 \ 

n jtk = exp(-0.7(j - fc)«7r)exp ( — 4000~^ ~ ^ ) ' 1 -^ k - m ( 75 ) 

with i = y/—l. As expected, the analytical curves match closely with the 
simulated curves. 

4. New Maximum Eigenvalue Distributions 

In this section, we shift attention to the distribution of the maximum 
eigenvalue of correlated non-central Wishart and gamma- Wishart random 
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Figure 2: Comparison of the analytical minimum eigenvalue c.d.f.s with simulated data 
points for correlated gamma- Wishart matrices with various dimensions and parameters. 
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matrices. As for the minimum eigenvalue distribution considered previously, 
once again the most direct approach of integrating the joint eigenvalue p.d.f. 
over a suitable multidimensional region seems intractable. To this end, we 
write the maximum eigenvalue A maa; (Y) of Y 6 as 

F max {x) = P (\ max (Y) < x) = P (Y < xl m ) (76) 

which allows one to deal purely with the distribution of Y, rather than the 
distribution of its eigenvalues. 

4-1- Correlated Non-Central Wishart Case 

For the non-central Wishart scenario, we deal with the matrix W with 
joint density given in fl3]). Thus, with ( ITS"]) , we have 



P (A ma:E (W) <x)= / w (W)dW 

'W<il m 



exp(— rf) 

f m (n)|S| n Jxi m -weni 

x F X (n; ©E _1 W) dW. (77) 



f |W| n - m etr (-S _1 W) 



Applying the change of variable W = xY with dW = x m2 dY in (1771) gives 

rf,mn ( _ ri \ flm 

P (X max (W) <x) = \ Pl V) / |Yr" m etr (-xS^Y) 
1 1 ; ; r m (n)|S|" J ' ' V ; 

x qFx (ra; x0S^Y) dY. (78) 

Expanding the hypergeometric function with its equivalent series expansion 
followed by using the reasoning which led to (1281) yields 



P (A_(W) < x) = C^tA y M [ " | Y|— etr (-x^Y) 

xtr k (aa H Y)dY (79) 
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where we have applied (cx H Yct) k = ti k (cxck^Y). This matrix integral seems 
intractable for arbitrary values m and n. In fact, this integral seems even 
more difficult to tackle than that which arises in the minimum eigenvalue 
formulation, i.e., Eq. ( 130]) . As the following theorem shows, however, we can 
obtain a solution for the case of 2 x 2 non-central Wishart matrices with 
arbitrary degrees of freedom. This is significant, because it presents the 
first tractable result for the maximum eigenvalue c.d.f. of correlated complex 
non-central Wishart matrices. 

Theorem 6. Let X ~ CN n ,2 (T, In ® S), where T G C nx2 has rank one, 
and W = X^X. Then the c.d.f. of X max (W) is given by 



x 2n exp(— rj) ^-v (xff) 



-ID \- yxfi) (fc) 
n !( n + l)! 2^( n ) fcife! 9 , - a :E- 1 ,«aH, n W (8Uj 

fe=0 



where 0_ s -i QCeH (0) calculated recursively via (TU)-(Tl\). 

Proof. Substituting m = 2 into ( 1791) . the proof follows upon application of 
Lemma |3j □ 

Remark 6. An alternative expression for can be obtained by employing 
the moment generating function based power series expansion approach given 



in fSd] . However, we have found that by employing that approach the final 
expression is more complicated, since it includes two infinite summations 
along with a recursive summation term. 

4-2. Correlated Gamma-Wishart Case 

We now turn consider the maximum eigenvalue distribution of gamma- 
Wishart random matrices. In this case, we deal with the matrix V with joint 
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density given in Thus, with ( 175|) . we have 

P (A maa; (V) <x) = K m , n x mn / |Yr m etr (-xQY) 

Jo 

x X F X (n - a; n; -arfSY) dY. (81) 

In the following theorem, we present a new exact closed form expression for 
the c.d.f. of the maximum eigenvalue of V for some particularizations of m, 
n and a. 

Theorem 7. Let V ~ rW 2 (n,a,S,n) with a > n > 2. Then the c.d.f. of 

A m ax(V) is given by 

2(a— n) min(fc,(o— n)) | 2 | 

^(x) = /C 2 ,„x 2 " ^ ^fci,i(^)^ fc (82) 

k=0 fc 1= [i] *=o 

ui/iere 

t-, / \ r 2 (2)r 2 (1/^^ + 2) ( £fc ,) 

^ fcl ' /(x) = V(u +A\ ^,S,^, 1 +2(°). ( 83 ) 

e^,! = 2/ci — k — 21, Uk lt i — n+ I + k — k\ — 2, k — (k\, k — ki) is a partition 
of k such that k\ G {0, 1, . . . , (a — n)} and |~|] < ki < min (k, (a — n)). The 
term (f'-xC^s,^ (+2(0) ^ s calculated recursively via < f7^) -< f77)j . 

Proof. Particularizing flHTT) tom = 2, a>n>2 and a G Z + and applying 
the zonal polynomial expansion (JSJ) gives 

F ma ,(x) = ^ 2 « 2( f; n) E [ - ( "; w)] " ( -f fc / l2 |Yr- 2 etr(-xQY) 
^ ^« [n]« k\ Jo 

x C K (SY)dY 
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Following the similar reasoning which led to (I56p . with some algebraic ma- 
nipulations we obtain (182]) . but with 



K K j(x) = [ 2 etr (-xQY) \Y\ Uk iHi £k ^(SY)dY. 
Jo 

This integrals is solved via Lemma [3] to yield ( I83p . □ 

Note that the c.d.f. result in Theorem [7] can be evaluated numerically for 
any value of n. Moreover, for specific values of n it often gives simplified 
solutions. Some examples are shown in the following corollaries. 

Corollary 4. Let V ~ TW 2 (n, n + 1, S, ft). Then the c.d.f. of X max (V) is 
given by 





n| 




n!(n + 1)! 




1 " 


S" 1 + n\ 


n+l 



X 



x ) = u , ■ ( i^i n + 2 ; -*Q) + r^Q,s,n(0) 



n 



+ (, + l l S )K + 2) ^'" +1; " + 3 -- lQ 0- <84) 
Corollary 5. Let V ~ TW 2 (n,n + 2,£,fJ). T/jen tfie c.d./. of X max (V) is 
given by 





n| 




n!(n + 1)! 


|S 


In 


E -1 + n| 


ra+2 



n 



x 2 (2 ) _ . 2|S|x 2 



2\S\x 3 (1) 



L F 1 (n + 2;n + 4;-xQ) . (85) 



(n+ l) 2 {n + 2) 
|S|V 

(n + l)(n + 2) 2 (n + 3) lJ 
Fig. [3] compares the analytical c.d.f. results for the maximum eigenvalue 
of gamma- Wishart matrices with simulated data, the matrix X and f2 are 
constructed as in ( 150]) and ( 1751) respectively. The analytical curves were 
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Figure 3: Comparison of the analytical maximum eigenvalue c.d.f.s with simulated data 
points for correlated gamma- Wishart matrices with various dimensions and parameters. 

computed based on Theorem [7J The agreement between the analysis and 
simulation is clearly evident. 

5. Conclusions 

We have derived new exact closed- form expressions for the c.d.f. of the 
extreme eigenvalues of correlated complex non-central Wishart and gamma- 
Wishart random matrices. We would like to conclude by emphasizing that 
these results provide the first tractable exact analytical results pertaining 
to the eigenvalue distributions of both complex non-central Wishart and 
gamma-Wishart random matrices with non-trivial correlation structures. Ob- 
taining tractable solutions for extreme eigenvalue densities for generalized pa- 
rameters (e.g., for arbitrary matrix dimensions) remains an important open 
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problem. 



Appendix A. Proof of Lemma [2] 



Proof. We start by factorizing x™ — xV^ and using x\ + 22 = tr(X) and X1X2 
I XI to obtain 



1 



n-2 
2 



tr(X)J](tr 2 (X)-4|X|cos 2 (^ 



Xi - x 2 



n-l 
2 



[ [ (tr 2 (X) -4|X|cos 2 



Next, recalling the generating function expansion 



i 


7Tj_ 




n 




)) 


n 





even n 



odd n. 



(A.l) 



(A.2) 



where denotes the ith. elementary symmetric function [34J of the parame- 
ters ip 2 , ■ ■ ■ 1 VVi}> an d using (IA.2[) in (IA. 1 [) along with some algebra, we 
obtain the result. □ 



Appendix B. Proof of Lemma [3] 

Eq. 3.23], we have 



Using 
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i 2 



x F t (a; a + 2; X) = K, / |Z| a - 2 etr (XZ) dZ 

'0 



(B.l) 



where K, = r2 ( a + 2 ) anc [ sft( a ) > 1. Following the proof of 28, Lemma 71, we 

l 2(0)1.2(4) V ' u J. E 

substitute X = A + By into ( IB. II) to yield 



iF x (a; a + 2; A + By) = /C / |Z| etr ((A + By) Z) cfZ (B.2) 

Jo 
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where y > 0. Expanding the term etr(BZy) gives 



iFi (a; a + 2; A + By) = /C V ^- / |Z| a - 2 etr (AZ) tr p (BZ) dZ. (B.3) 

Now, we aim to establish a power series expansion for \Fi (a; a + 2; A + By) 
around y = 0. To this end, denote 



0A,B,a(y) = i-Fi (a; a + 2; A + By) . 



(B.4) 



We then have 



L F X (a; a + 2; A + By) = £ ^ B , a (0). 



p=0 



(B.5) 



Equating the coefficient of y p on both sides of flB.3j) and (IB .5|) gives (jUJ). 



Following 
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45], we can express the confluent hypergeometric function 



(B.6) 



of a matrix argument in the determinant form 

, / s A A ,B,a(j/) 

"-a,b (y) 

where h^p{y) = X\{y)—X2{y)- Since we are interested in obtaining AB a (0), 
we may rearrange (IB.6P such that 

^A,B(l/)0A,B,a(y) = A A ,B,a(y) 



and apply Leibniz's rule 19] for the fcth derivative of a product to obtain 



E = Ai fc) Bja (y). 

i=o ^ 



(B.7) 



After rearrangement of terms we obtain the following recursive formula 

C,«W = r — 73 ( R8 ) 
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which, upon evaluating at y = 0, gives (JTUJ). 

What remains is to evaluate the successive derivatives fe^ B (0); equiva- 
lently, x^\o) and x% (0). To this end, we use the relations 

zi(y) + ar 2 (l/) = tr (A) + ytr (B) 

(B.9) 

xi(y)x 2 (y) = |A + By| = |A| + |A|tr (BA x ) y + |B|y 2 . 
Evaluating the first derivative of (1B.9j) with respect to y at y = gives 

a;S 1) (0) + 4 1) (0) =tr(B) 

(B.10) 

x 2 (0)x l i\0) + xi(0)4 1} (0) = |A|tr (BA" 1 ) 

which upon solving for x{ (0) and 4 (0) gi yes the corresponding results in 
(Q3J) and (TT5l) (i.e., j = 1). Taking the second derivative of (1B.9j) . followed by 
similar calculations as before, gives the case corresponding to j = 2 in (|T%|) 
and (|15p . The remaining case, j > 3, is more challenging. To proceed, let us 
take the jth derivative of (IB.9j) for j > 3 to obtain 

x?' ) ( 3 /) + xJ' ) ( 3 /) = 

(B.H) 



fc=0 

where we have again used the Leibniz's formula to obtain the jth derivative 
of the product X\(y)x2(y)- After some rearrangement of terms followed by 
evaluating the resultant derivatives at y — gives 

x[ j) (0) + x { 2 j \o) = 

/n\ (B.12) 
xS j) (0)x 2 (0) +4 j) (0)x 1 (0) = -J2 ( l)^~ fc) (0)4 fc) (0)- 

fc=i ^ ' 

These simultaneous equations can easily be solved for Xi\o) and x 2 \o) to 
yield the results in fTHj) and f|T5|) . □ 
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Appendix C. Proof of Lemma SI 



For Z G and R G T-L m with rank one, we have from [37], Eq. 6.1.20] 

f etr (-ZX) |X| a ~ m C r (XR) dX = (a) t f m (a)|Z|~ a a (RZ^ 1 ) (C.l) 

where 3?(a) > m — 1 and r is a partition of £. Following the proof of 28 . 
Lemma 7], let us now select Z such that Z = A + yl m , where A G and 
y > 0. Substituting this specific value of Z into ( IC.lj) and choosing R such 
that R = rr^ where r G C mxl , yields 

/ etr (—AX - yX) |X| Q ~ m tr* (XR) dX = 

(a) t f m (a) | A + yl m \~ a (r 11 (A + ylj" 1 r)\ (C.2) 
Moreover, we can expand the term etr (— yX) to obtain 

£^T"/ /M-AX)tr fc (x)|xr m tr<(XR)dx = £(?/) (C-3) 

where 

£(</) := (a) t f m (a) | A + yLnp (r H (A + yU" 1 r)' . (C.4) 

Now, we seek a power series expansion for the real- valued function £(?/) 
around ?/ = 0. Equating the coefficient of y with that on the left-hand 



side of ( 1C.3j) will then give the desired expression. 

We require £^(0). To evaluate this, we start with the eigen-decomposition 
A = USU^, where U G C mxm is unitary and £ = diag (at, a 2 , . . . , <f m ), to 
obtain 

\hi 



> t 

12 



{(») =(a)(T m (a) x '" * (C.5) 



rife+ 
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and U^r =: h = (hi h 2 . . . h m ) T . It is not difficult to see that we can 
have a convergent power series if we select y < min (oS, <t 2 , . . . , a m ). Finally 
equating £^(0) with the coefficient of y on the left-hand side of flC.3j) with 
J2 ^ = tr ( r H (A" 1 )" r) = tr (R (A" 1 )") gives (©. □ 



i=l 



Appendix D. Proof of Lemma [5] 



We combine (1C.3j) and (1C5[) for the case m = 2 and apply the relation 



cr 2 |/ii| +0i|/i 2 | — |A|tr(RA ) to arrive at 



p=0 

where 



y>(_JJV /" etr(-AX)tr p (X)|X| a - 2 tr < (XR)rfX = C(y) (D.l) 



{ - M =7>- ^ ( D . 2) 



with p = W.f.taMR) ta^ ffld & _ |A|tr (RA">) 



|A|*+ a ' " 2^/\A\ tr(R) 
Our objective is to obtain a power series expansion for ((y) around y = 0. 

To this end, we may use the generating function definition of ultraspherical 

polynomials^ to write 



^ +?f . + r = fsm. (D ,) 

|A| ,/lAI ) f l A l» 



Now we may use (ID.3I) in ( ID. 21) with binomial theorem to obtain 



c(i/) = ^EE (O^^tai/ 1 ^ • (D - 4) 

n=0 1=0 ^ ' I I 2 



2. 



3 Ultraspherical polynomials can be denned through the generating function as [2|, Eq. 
6.4.10] (l-2xr + r 2 )- x = J2n=o C n(x> n - 
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Since the desired general form of the expansion is 

C(y) = E L -r^' (D - 5) 

what is left is to determine the coefficient A p using (ID.4|) . To this end, we 
must collect the coefficients of y p together. Since (1D.4|) contains a finite 
inner summation, we have to consider two cases depending on the value of 
t; namely p < t and p > t. When p <t, the summation indices are selected 
from the set l,n — {0, 1,2, ... ,p} such that I + n = p. In the case of p > t, 
the summation indices are selected from the sets I = {0,1,2, ... ,t} and 
n — {p — t,p — t + l,p — t + 2, . . . ,p) such that I + n = p. Putting these 
together, we come up with a new set 

k = {0,1,2,..., min(p, t)}, n = p-k (D.6) 

which embraces both cases. Thus, the coefficient A v can be written as 

A P =v P \ y: ^)ilY k -7vM ^ 

k=o W l A l 2 

where we have used the fact that z) = (— l) n C^(z). Using this, equating 
the coefficient of y p in (ID. 51) and (ID. If) concludes the proof. □ 

Appendix E. Proof of Lemma [6] 

Before proceeding, it is worth mentioning the following relation 

C 1|1)0 (X) = |X|tr (X- 1 ) = |X|C 1A0 (X- 1 ) (E.l) 



where X e Wl^ . Also, for Z e T-L£ , we have from 



46l Eq. 3.10] 



etr (-ZX) |X|C li0)0 (X- i )dX = r 3 (4)|Z|- 4 C 1A o(Z). (E.2) 

xeni 
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Following the proof of [28|, Lemma 7], let us substitute Z = A + Ry, for 
y > 0, into (IE. 2 j) to obtain 



I etr (—AX - RXy) iXld^X-^dX 



r 3 (4)|A|- 4 \l 3 + A- L Ry\ tr ( A + Ry) 



(E.3) 



Since R is unit rank, A X R is also unit rank, and therefore (1E.3|) can be 
written as 



f 



etr (—AX - RXy) iXlCVo^X^dX 

= f 3 (4)|A|- 4 tr (A + Ry) 1 F (4; -tr (A^R) y) (E.4) 



where we have used the relation 1/(1 + z) n = iF (n; —z). Now, since y is 
arbitrary, we select y < l/tr(A _1 R) to obtain a power series expansion for 
the right-hand side of (IE. 41) as 



etr (—AX - RXy) iXlCi^X^dX 



f 3 (4)|A|- 4 tr(A+R 2/ )J]^tr' (A^R) (-yf. (E.5) 



t=o 



t! 



Finally, expanding the left side of (1E.5|) as a power series of y followed by 
equating the coefficient of (— y)* on both sides with some manipulations con- 
clude the proof. □ 



Appendix F. Proof of Lemma [7] 



We first solve 



/ etr (—AX) tr p (BX) |X| a ~ 2 C T (X)dX 
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Subsequent application of the basic property ^ T C T (X) = tr*(X) will then 
yield the desired result. 



Let us begin with the following matrix integral 



33, Eq. 6.1.20 



/ etr(-ZX) |X| a - 2 CV(X)dX = f 2 (a)[a] T |Z|- a C T (Z- 1 ) (F.l) 

where ZeH 2 + and U{a) > 1. Selecting Z = A + By, where A,B e Ut, 
fIF.ip becomes 

f etr(-AX-BXy)|X| a ~ 2 C r (X)rfX = C(y) (F.2) 

where 

((y) = f 2 (o)[o] T |A + By\- a C T ((A + By)' 1 ). (F.3) 



Since the left-hand side of (1F.2|) can be expanded as a power series in y, the 



remaining task is to find a power series expansion for the right-hand side of 
( lF.2p . i.e., C(y), so that the coefficient of y p can be compared on both sides. 
To this end, we expand the zonal polynomials in flF.3j) using [37, Eq. 6.1.12] 
to obtain 

C(y) = f 2 (a)[a] T ^|^±^|A + By\-^h (FA) 

*l-*2+!_ *l-*2+ 1 

where 7 = — an d are the eigenvalues of A + By. At this 

point, observe that since (ti,£ 2 ) is a partition of t, we can write t 2 = t — t\, 
where [|] < ti < i. With this observation and the aid of Lemma [2], we then 
obtain 

7 = (-l)*4\ T tr £t ^ (A + By) |A + By\- {a+£t ^ . (F.5) 

i=0 
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Next, with the binomial expansion we get 

r2t 1 -t-i"| 
I a I e H,i 



C(y) = K tl f 2 (a) E(-l) 4 4 i e[(^' t )tr- 1 - fe (A)tr fc (B) 

A + By\- {a+£t ^y k . (F.6) 



i=0 fc=0 

X 



where •= t' ( a )«i("-i)*i-i(2*i-*+i) 

We now aim to obtain a power series expansion for |A + Byp < - a+e * 1 ' ) in 
terms of y. To this end, we may express 

|_AJ-(a+eti) 



|A + By| 



-(a+£ tl ) 



(l + 2/3v^B|y + |A-iB|/) a+£tl 

oo 

| A |-(a+ etl ) ^ ia^BIS C +£tl (/?) (-y)" (F.7) 



n=0 

tr(A- 1 B , l 



where /3 := — , with \y\ < , 1 Here, to obtain the last equality in 

2^/|A" 1 B| m -y/lA-iBI H J 

( 1F.7|) . we have exploited the generating function definition for ultraspherical 
polynomials |2.!, Eq. 6.4.10]. 

Incorporating ( IF. 71) into ( 1F.6I) gives 



((y) = K tl T 2 (a) J2 E(-l) t+fc 4 i e[(^Vr- 1 ^ fc (A)tr fe (B) 



i=Q fc=0 

oo 



X 



| A |-(«+^) lA^Blf C^ +etl (/3) (-y) n+fc . (F.8) 



n=0 



Since we are interested in the coefficient of (—y) p , we have to re-sum the 
above series to collect all terms having power {—y) p . A careful inspection of 
the above equation reveals that we can select n = p — k and the upper limit 
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of k as mm(p,£t lt i). Thus, we have after some manipulations 

r2t 1 -i-l"| 
OO | 2 I 

ay) = K tl t2(a)\A\- a J2 £ B ^A-y) p - (F.9) 

p=0 i=0 

Now, equating the coefficient of {—y) p with the corresponding coefficient in 
(jK2j) yields 

r 2ij^ — i — 1 "1 

/. I 2 | 

/ etr(-AX)tr p (BX)|X| a - 2 a(X)rfX = Z ilP !f 2 (a)|Ar a V B Tpi . 

Finally, using the basic property ^ T C r (X) = tr*(X) along with the fact that 
Y^ T = Z/t 1= r*l gi ves the desired result. □ 
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